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COLLISIONLESS DAMPING OF FAST MHD WAVES IN MAGNETO-ROTATIONAL WINDS 
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ABSTRACT 

We propose collisionless damping of fast MHD waves as an important mechanism for the heating and ac- 
celeration of winds from rotating stars. Stellar rotation causes magnetic field lines anchored at the surface to 
form a spiral pattern and magneto-rotational winds can be driven. If the structure is a magnetically dominated, 
fast MHD waves generated at the surface can propagate almost radially outward and cross the field lines. The 
propagating waves undergo collisionless damping owing to interactions with particles surfing on magnetic mir- 
rors that are formed by the waves themselves. The damping is especially effective where the angle between 
the wave propagation and the field lines becomes moderately large (~ 20 to 80°). The angle tends naturally to 
increase into this range because the field in magneto-rotational winds develops an increasingly large azimuthal 
component. The dissipation of the wave energy produces heating and acceleration of the outflow. We show 
using specified wind structures that this damping process can be important in both solar-type stars and massive 
stars that have moderately large rotation rates. This mechanism can play a role in coronae of young solar-type 
stars which are rapidly rotating and show X-ray luminosities much larger than the sun. The mechanism could 
also be important for producing the extended X-ray emitting regions inferred to exist in massive stars of spectral 
type middle B and later. 

Subject headings: magnetic fields - plasma - magnetohydrodynamics - stars : winds - waves 



1. INTRODUCTION 

A wide range of stellar classes are thought to have hot plas- 
mas with temperatures greater than 10 6 K in their atmosphere. 
Main sequence stars with low to intermediate mass, including 
the sun, are known to possess hot coronae, and stellar winds 
emanate from coronae by thermal expansion most likely en- 
hanced by magnetic wave energy deposition. Magnetic fields 
generated by dynamo mechanisms in these stars are likely 
brought to the surface by magnetic buoyancy in the convec- 
tion zones JParker Hi 9660 . and the fields play a major role in 
the structure of the outer atmosphere. In a static atmosphere, 
the gas pressure, p, falls exponentially with a scale height H, 
but the magnetic field, B, generally decreases more slowly, as 
a power-law with radius r, to satisfy conservation of magnetic 
flux. As a result, the magnetic pressure, B 2 /8w, can eventu- 
ally dominate thermal pressure from a certain radius outward, 
even though the thermal pressure dominates in regions nearer 
the photosphere. Thus in the outer regions the dissipation 
of the magnetic energy can greatly influence the energetics 
and dynamics of the atmosphere and wind, producing in some 
way, the heating of coronae and the wave acceleration of stel- 
lar winds. We intend to explore a wave damping process that 
can produce the heating and acceleration. 

Stellar rotation is also closely related to the coronal activ- 
ities. The observed X-ray luminosity of low to intermediate 
mass stars increases monotonically with decreasing rotation 
period, P for P > several days, and th en is saturated for 
smaller values of P JPizzolato et al.l2003l) . It is likely that the 
stronger differential rotation in the interior leads to the gener- 
ation of a surface magnetic field that is stronger, and a more 
active corona is formed. However, even qualitative arguments 
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are still lacking as to how the heating and acceleration by the 
wave damping actually occurs in the atmospheres of rotating 
stars. 

It is now becoming more evident that the heating and driv- 
ing effects of magnetic fields is not confined only to cool stars, 
but can be important in the more massive early type stars as 
well. For many years it had been believed that dynamo gen- 
erated magnetic fields would not exist on hot stars, and that 
the winds from these objects were produced only as a result 
of radiation pressure gradients. Unlike the cooler stars, the 
hot massive stars do not have an outer convection zone, so the 
mechanism for the rise of magnetic fields and the source of 
mechanical energy for coronal heating are not present. How- 
ever, early UV spectra of the stars showed lines from a noma- 
lously high ionization stages dLamers & Mortonl 19761) . This 
superionization could be expla ined as resulting from X- ray 
ionization by the Auger Effect (Cassinelli & 01s onlfl979i) . It 
was soon realized from satellite observations that essentially 
all O and B stars are X-ray sources (Seward et al. 1979; 
Berghofer et al. 1997), bu t the nature of the source o f the 
X-rays remained a mystery. ICassinelli & Swankl i 19831) pro- 
vided arguments that perhaps there are both shocks in the 
winds (as proposed by Lucy 1982), as well as hot coronal 
regions at the base of the winds. This idea that there are two 
contributors to the X-ray emission has persisted (Waldron and 
Cassinelli, 2000). Also, in the case of B near main sequence 
stars, Cohen et al.(1997) found that the X-ray emission mea- 
sure was comparable-to or larger than that of a theoretically 
predicted wind, and hence the coronal component of the X-ray 
emission could be dominant for these stars. As for the mag- 
netic fields in the outer atmospheres and winds, some direct 
measurements have been made in near main sequence B stars 
with fields of about 300 Gauss (Donati et al. 2001). Char- 
bonneau and MacGregor (2001) have shown that fields can be 
generated by a dynamo operating at the interface between the 
convective core and the radiative envelope in early type stars. 
MacGregor and Cassinelli (2003) have shown that such fields 
can buoyantly rise through the radiative zone to the surface of 
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hot stars. In addition there is a wide variety of observational 
evidence that hot stars have magnetic fields (Henrichs 2001). 

Thus, there appears to be a greater similarity between hot 
and cool stars than had been envisioned, and it is both timely 
and useful to consider the processes by which magnetic fields 
can heat and accelerate winds across the HR diagram. Here 
the focus is on structures in which the field in the wind is 
configured in the form of a spiral, as would naturally arise 
from the rotati on of a star with a mag netic field rooted in its 
surfac e (fig.Q) JWeber & Davis H967llBelcher & MacGregorl 

MM- 

Magnetohydrodynamic (MHD) waves are widely regarded 
as promising mechanisms that play a role in the heating and 
acceleration of the outer atmospheric plasmas. Surface turbu- 
lence and transient activities (e.g. Sturrock 1999) excite MHD 
waves at the surface of the sun. It is natural to expect that 
similar excitation phenomena occur on other stars that have 
surface magnetic fields. Outwardly propagating MHD waves 
carry Poynting flux energy, and any dissipation of this flux 
directly leads to conversion from magnetic energy to thermal 
and kinetic energy in upper atmosphere. There are numer- 
ous studies of wave heating in the solar corona / wind (e.g. 
Belcher 1971) and in stellar coronae / winds (e.g. Hartmann 
& MacGregor 1980). M HD waves are d iscussed in regards to 
pulsar winds as well fLyubarskv 120031) . 

A flow of wave-like motions created by turbulence in a lo- 
calized source can be decomposed into three types of MHD 
waves: fast, Alfven , and slow modes (see Cho & Lazarian 
2003). In a magnetically dominated plasma, the slow mode 
essentially corresponds to an acoustic wave, hence, it prob- 
ably does not contribute the dominant amount of heating to 
the plasma in the abse nce of other strong driving mechanisms 
(as discussed in i|5.3> . The fast and Alfven modes, however, 
can play a dominant role since their wave energies can easily 
exceed the thermal energy. These two modes propagate quite 
differently. The Alfven waves propagate as transverse waves 
and only along the field lines, while fast waves can propa- 
gate almost isotropically. The properties of fast mode wave 
do depend on the direction of the wave relative to the mag- 
netic field line and thus can change with radial distance from 
the star. The component of the fast mode wave which is trav- 
eling along the field line is transverse, as is the case for the 
Alfven mode, but the fast wave becomes increasingly com- 
pressive as the angle, 6, between the wave propagation vector 
and the field line increases. 

Some of the properties of the wave propagation are affected 
by the rotation of the star because the rotation causes the mag- 
netic field lines that are rooted on the star to form a spiral 
pattern as described above. In the solar wind, the fiel d is usu- 
ally p ictured as developing a so-called Parker spiral iParker I 
119631) . This spiral arises in spite of the relatively slow rotation 
of the Sun. More tightly wrapped spiral patterns should form 
in the winds of stars that rotate faster, such as younger solar- 
mass stars and massive stars. Thus, it is important and inter- 
esting to study outward propagation of the fast MHD waves in 
magneto-rotational winds of rapidly rotating stars in general. 

Let us consider outgoing fast waves emergent from the stel- 
lar surface. The fast waves will travel almost radially with- 
out being refracted provided that the magnetic energy in the 
plasma dominates the thermal energy. However, as the wave 
travels away from the surface, the angle of propagation, 6, in- 
creases monotonically (fig. Q and the character of the wave 
changes accordingly. The way in which the dissipation of the 
waves is influenced by this variation of 6 is described in (SJ3, 




FIG. 1. — Schematic picture of fast MHD waves propagating outward in 
magneto-rotational wind, illustrating the angle 9 between the propagation 
direction and the field at one point in the outer atmosphere. 



where we also discuss how the stellar rotation plays a domi- 
nant role in controlling the wave dissipation and subsequent 
heating of the stellar atmosphere (f|4}. 

An important aspect in the discussion of stellar plasmas 
is that they are generally collisionless except in regions very 
near the surface. This is the case because the mean free path 
for a Coulomb collision, 

(T/10 6 K) 2 

Up = 9.38 x 10 7 cm-^-g (1) 

(n/10 8 cm A ) 

is larger than wave length. Thus it is essential for us to con- 
sider collisionless processes in the wave dis sipation. The fast 
mode waves suffer collisionless damping (lGinzburg|ll961t 
lBarneslH%6T) . which is discussed in the context of the heat- 
ing of the solar wind plasma by IB arnesl J 1 9681 1 1 9691) . Inter- 
estingly, the damping rate sensitively depends on 9 (S0. To 
date, however, this damping process has not been applied to 
magnetically dominated winds from rapidly rotating stars. It 
is the goal of this paper to investigate the importance of col- 
lisionless damping for both the heating and acceleration of 
stellar coronae in which the field has the configuration of a 
magneto-rotational wind. 

2. DAMPING OF FAST MHD WAVES 

In this section, we summarize the physical process of col- 
lisionless damping of the fast MHD waves, and we introduce 
the damping rate. It is useful to compare the collisionless 
damping proces s with the one that has often be en considered 
in the literature dHollweg 111981 ISuzukill200l : the damping 
of waves through the formation of a fast shock train by the 
steepening of the wave fronts. We find circumstances in coro- 
nae under which the collisionless damping can become more 
important than the shock steepening process. 

As discussed in Ginzburg (1961), the nature of collision- 
less damping of a wave can be considered analogous to the 
creation of radiation by charged particles in magnetic field. 
The charged particles can emit electromagnetic waves both 
by their acceleration (which produces cyclotron radiation) and 
by the Cherenkov effect. By the inverse of these two pro- 
cesses charged particles can absorb the radiation. Similarly, 
plasmas can absorb energy from plasma waves by a gyrores- 
onance process or by a process analogous to the Cherenkov 
effect. It is the latter that is of interest to us here. It leads to 
a dissipation of the plasma waves by a wave-particle interac- 
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Fig. 2. — Propagation of fast MHD waves for two different angles, 9. 
In a low (3 plasma, the displacement vector, £ j , of the fast mode is almost 
perpe ndicular to underlying magnetic field, Bq t Cho, Lazari an. & Vishniad 
2003). This illustrates that the distance between the two oscillating field lines 
(solid lines) changes more in larger 9 case than in the smaller 9 case. This is 
the case even though, as illustrated, the absolute displacements, \£f\, are the 
same. Particularly, the closest distance between the two field lines (indicated 
by dotted both-side-arrows) is smaller in the large 9 case, leading to a more 
reflective mirror. Hence, a larger number of particles are locked into the 
surfing on the mirrors. This in turn leads to a higher rate of wave dissipation 
by particle-field interactions. 



tion, which results in heating and acceleration of the particles 
that compose the plasma. 

The resonance of a wave with the gyro-resonance frequency 
of ambient thermal ions causes a damping of those wave 
modes that have a frequency close to the ion-cyclotron fre- 
quency (e.g. Leamon et al.1998). Thermal particles can also 
be accelerated by an oscillating parallel electric field, as is the 
case for Landau damping of a wave. A similar mechanism 
called transit time damping (TTD) opera tes if the electr ic field 
is replaced by moving magnetic mirrors, barnesl Jl966') intro- 
duced this idea to expl ain the heating of the solar wind, and 
l lYan" & Lazarian 12004ft applied to studies on cosmic ray scat- 
tering. In this paper, we consider the TTD of fast MHD waves 
in rapidly rotating stars. We can reasanably assume a magnet- 
ically dominated medium, commonly called a low (3 plasma. 
As usual, (3 is the the ratio of the gas pressure p g = pa 2 , to 
the magnetic pressure, pi, = B 2 /8ir, 



Pb 



8-irpa 2 
~B 2 ~ 



(2) 



where B is the magnetic field strength, p is the ambient gas 
density, and a is the isothermal sound speed. 

The TTD is caused by interactions between magnetic mir- 
rors associated with the fast waves (Figure |2ji and particles 
that have a parallel velocity, v\\, similar to the pattern speed of 
the mirrors, that is, 

v\\ ^ u/k\\ . (3) 

This is the Cherenkov condition, and u) and fen are the wave 
frequency and the wave number parallel to the magnetic field 
lines. Particles which satisfy this condition surf on the mirrors 
and thus interact with the waves. The damping rate 7^ of fast 
waves which propagate with an angle 6>(~ 1 radian) relative 
to the field lines in a plasma of a given /3(<C 1), and for the 



Maxwell-Boltzmann particle distribution is dGinzburg 119611) 

sin 2 9 
4 cosv 



m e m e 1 

ex P( a Ya> + 5cx P(~ o Ta> 

m H tuhP cos 2 9 p cos 2 9 



,(4) 



where mjj and m e are the hydrogen and electron masses, and 
this expression is strictly true only for linear waves in which 
the velocity amplitude is smaller than the phase speed. Note 
that 9 ch anges owing to the turbulent fluctuations in actual sit - 
uations (Lazarian & Vishniac lfT999l lYan & Lazarian 112004 . 
However, since we are considering small relative amplitude, 
SB/B, and are primarily interested in the average properties, 
we choose here to treat 9 as the average angle between the un- 
derlying magnetic field and the wave propagation direction. 

The Upper panel of fig. [3] presents 7<j divided by uj as a 
function of the angle 9 for two values of (3; (3 = 0.1 and 
0.01. This collisionless damping is shown in comparison with 
the damping rate associated with a shock train that is usually 
used for such discussions, (see the Appendix for a discussion 
of the the damping of waves by the shocks). The rate of the 
collisionless damping is independent of wave amplitude, SB, 
because it is a linear process, while it is larger in plasma which 
have a (3 value not too far below unity. This is because the 
parallel speeds of the surfing particles needs to be comparable 
to fast mode speed for the resonance to occur. If (3 is too low 
(i.e. the sound speed is too low), only an insignificant fraction 
of the thermal particles will have the large values required for 
the condition given in eq.(3). 

Collisionless damping is also quite sensitive to the angle 
9. The damping rate increases more rapidly than the steepen- 
ing in the range, 9 < 80°, but then it suddenly decreases as 
9 comes closer to 90°. This dependence on 9 can be under- 
stood by considering two competing factors. In general, the 
damping increases with 9 because a magnetic mirror tends to 
be formed in a more efficient way, as is illustrated in Figure 
13 A small separation of the field lines is an indication that 
a magnetic mirror has formed. The closest distance between 
the field lines become smaller for larger 9 so that the mir- 
rors become more reflective even though the amplitudes of 
the fast waves are the same. Thus, more particles are inter- 
acting with the wave to cause the dissipation. However, when 
9 approaches to 90°, fen becomes small, and the resonance 
condition (eq. 3) requires a large v\\ . Hence, most thermal 
particles will no longer be in resonance with the wave, so the 
damping rate decreases sharply in the 9 near 90° regime. 3 

On the other hand, the damping caused by the shocks does 
depend on SB because it is a nonlinear process, and is almo st 
independent of the plasma (3 value (see Appendix eq. JA^H . 
So, we adopt (3 = 0.01 for the two cases of the shock damp- 
ing. The damping rate by the shock steepening also increases 
with 9 as in the collisionless process, because fast waves con- 
tain a longitudinal, i.e. compressional, component as 9 in- 
creases, and they become more dissipative. However, the 
variation is more weakly dependent on 9 than is the case of 
collisionless damping. 

The conditions under which the collisionless process ex- 
ceeds the shock steepening process is further illustrated in the 



The former effect appears in the 



term in eq.j^J which increases 



with 9, while the latter effect appears in the exponential terms in eq.JJJ which 
both approach as 9 90°. 
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Hnlfclnnle^, ft = 0' 1 

Collisionless, /S=0.01 

Steepening, 6B/B = 0.1 

Steepening, 6B/B=0.01 




FIG. 3. — upper: Comparison of damping rate of the collisionless process 
with the steepening as a function of 9. lower : Regions in 9 — j3 where 
the collisionless process dominates the steepening for dB/B = 0.1 (cross 
hatched) and dB/B = 0.01 (single hatched). 



lower panel of fig. [3] which shows on a 9 versus (3 plane for 
two cases, SB/B = 0.1 (cross hatched) and SB/B = 0.01 
(single hatched). Note that the collisionless process tends to 
become more important for larger values of 9, except when 9 
is near 90° and/or when the (3 value is too small. As is shown, 
a fast wave with amplitude 6B/B = 0.01 will dissipate pre- 
dominantly by collisionless damping over a wide range of 9. 
Even for the case in which 5B/B = 0.1, the collisionless 
damping is seen to be more important for (3 ~ 0.1 in an an- 
gular region of 60° < 9 < 90°. These results show that there 
are interesting physical conditions under which collisionless 
damping should, from now on, be taken to be the primary pro- 
cess for wave dissipation. 

These results have shown that the collisionless damping can 
be comparable-to, or more important than, shock steepening 
for the damping of fast waves. In fact, collisionless damping 
could dominate even more-so. This is because shock steepen- 
ing could be smaller than we have assumed, for the following 
reasons. First, when a sinusoidal wave is excited initially, it 
takes some time for it to steepen to form the shocks before 
the shock dissipation process can begin. Whereas the col- 
lisionless damping works immediately, even for the case of 
sinusoidal waves. Second, the collisionless effects can them- 
selves modify the shock dissipation. In general astrophysical 
circumstances, each shock in a shock train is considered col- 
lisionless, so the plasma is not thermalized by Coulomb colli- 
sions, but rather the plasma is randomized because of fluctu- 
ations in the magnetic field. In such situations, particles with 
energy larger than average can move across the front more or 
less transparently and the conversion from kinetic energy to 
heat would be less effective, and hence also the shock dissi- 
pation would be less. We conclude from this section that our 
results indicate that it is important to account for collisionless 
damping of the fast modes for realistic and interesting astro- 
physical situations. 

3. MODEL FOR MAGNETO-ROTATIONAL WINDS 



Let us consider the outward propagation of the fast MHD 
waves for the case of magneto-rotational stellar winds. Fig- 
ure^depicts propagation of the fast MHD waves in the rota- 
tional winds seen from the pole. Here we adopt a simple split 
monopole configuration for the magnetic fields and assume 
that the magnetic axis coincides with the stellar spin axis. Spi- 
ral magnetic fi elds are formed in the equatorial region by the 
stella r rotation JWeber & Davisll967tlBelcher & MacGregorl 
We can also reasonably assume that the ideal MHD 
condition holds in the stellar atmosphere. These are common 
assumptions for the equatorial region of a magneto-rotational 
wind. 

The plasma itself cannot move across the field lines. How- 
ever, the fast MHD waves can propagate almost isotropically 
in a low (3 plasma. In particular, they can propagate radially, 
crossing the field lines as illustrated in the figure. Also shown 
is the increase in the angle 9 as the wave travels away from 
the surface. Where the angle 9 is large, the fast waves suffer 
collisionless damping, as we discussed in the previous sec- 
tion as long as the wave length, A, is short enough to satisfy 
collisionless condition, 



A < I 



mfp ■ 



(5) 



For example, if we consider a fast mode wave with period of 
Is, with a phase speed is 1000 km/s traveling in plasma with 
a temperature of 2 x 10 6 K, the collisionless condition that 
must be satisfied for the particle density is n < 10 s cm -3 . 
In the case of the sun, this density condition holds even in 
the inner corona. As a result, the wave dissipation process 
directly leads to heating and acceleration of wind plasma and 
may have an effect on the wind structure. 

In this section, we present a model to explore the colli- 
sionless damping of the fast mode for the case of magneto- 
rotational winds. As we are just beginning to develop an un- 
derstanding of the basic collisionless damping process, it is 
too soon to be aiming for a complete and self consistent pic- 
ture. For example, there are other processes such as shock 
dissipation argued in the previous section that may still be 
important. We choose to estimate the energy and momentum 
that is transfered by the fast wave dissipation in the case of 
a pre-specified background wind. Then we can discuss the 
likely influence on the wind energetics and dynamics. 

3.1. Approximated Wind Structure 

Here we describe a practical method for constructing 
approximate wind structures in the equatorial plane of 
the magneto-rotational winds, as discussed in Lamers and 
Cassinelli (1999). We assume that the wind is isothermal at a 
temperature T, and we assume that all the physical quantities 
depend spatially only on the radial distance, r in the equatorial 
plane. For the split monopole magnetic field, the conservation 
of magnetic flux, V • B = 0, fixes radial component of the 
magnetic field, 

B r = H^J B rfi , (6) 

where R+ is stellar radius and B r ,o is radial field strength at 
the surface where r = i?*. It is empirically known that radial 
velocity distribution of stellar winds can be represented by the 
"beta velocity law" (here we use the power rj to avoid confu- 
sion with the plasma (3 ratio that we have used throughout), 
thus 

I 1 ) , (7) 
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where r\ ~ 0.5-3. If the rotation is slow, terminal velocity, 
Vrx,, is of order the escape velocity, u CS c,o — *J 2G r 1 * » where 
M* is stellar mass. 

In rapid rotating stars, w GSCi o may not a good indicator of the 
terminal velocity, as Voq > v O sc,0 because the wind is signif- 
icantly affected by acceleration from the magneto-rotational 
forces. For these cases, Lamers and Cassinelli (199 9) show 
JhatVoc is better estimated by the Michel velocity (tMichell 
EH), which is defined as 



Mr 



1/3 



Mr 



1/3 



(8) 



where O is angular speed of stellar rotation and Me is equa- 
torial mass flux rate multiplied by Anr 2 and we use, 



M E = 

In this paper we adopt 



4npv r r 2 (= const.) 



Voo = max(u esC! o, «m)- 



(9) 



(10) 



It should be noted that Voo can be estimated by the surface 
values since both u OS c,o an d v m are determined at r = R±. 
From here on, we follow convention and call cases with 
"osc.o > "Mi w OSCj0 ~ v M , and v cscfi < v M , the "slow", 
"moderate", and "fast" magneto-rotators, respectively. 

We can determine azimuthal field component, B§, for given 
O by the induction equation as 

V *~ rQ . (11) 



D,f, 

B r 



Azimuthal velocity, v^, is derived from conservation of angu- 
lar momentum as 



n 



y A,r 



V<f, 



where C is specific angular momentum, v\, r 



u A,r 



(12) 



/4irp 



dial Alfven speed, and ta is Alfven point at which v r — i>A,r- 
Thus, we have set the structures of the magneto-rotational 
winds for given M*, R+, Me, -B r .o an d So, where Sq is spin 
normalized by Kepler velocity, fix = 



GJJ, 



9 - " 



(13) 



3.2. Wave Propagation 

We introduce formulation describing propagation of the fast 
MHD waves with dissipation in the magneto-rotational winds 
under the WKB approximation. We only consider magnet- 
ically dominated i.e., low (3 plasmas. Accurately speaking, 
refraction affects the wave propagation; the fast wave prop- 
agates toward a region with smaller phase speed, V{ , defined 



as 



B 



vt 



1 



Air pa 2 
~B 2 ~ 



(1 



Am par 
~B 2 ~ 



AB 2 Anpa 2 
~B 2 W 



1/2 



(14) 



4 Although adiabatic sound speed should appear here, we use isothermal 
sound speed, a, as we assume the isothermal atmosphere. 



which corresponds to a direction away from the field line, 
and where B = < B 2 + B 2 . However, in the case of a 



low P plasma we can reasonably assume that the fast waves 
propagate radially without suffering the refraction, because 



Vi 



V47TP 



is almost independent of 9. 



The energy flux of MHD w aves i n fluids moving with ve- 
locity, v is written as dJacauesI 19771) 



V g £ w + v-p y 



(15) 



where £ w is wave energy density, V g is group velocity, and 
p w is wave pressure. For a wave amplitude, 8v, the energy 
density is £ w = \p8v 2 . The energy flux, Fw, changes ac- 
cording to work done by p w and damping energy loss rate, 
— Q w , so that in steady state, 



V • F w -v(V Pw ) = -pQ v 



(16) 



Note that Q w works as a heating term for the surrounding 
plasma. We are focusing on fast MHD waves in the low (3 
plasma. Thus the radial components of the pressure gradient, 
VjPfw i an d m e energy flux, i*f w , can be expressed as dJacauesI 
1977) 



(V.p fw ) r = (l+sm 2 6) dl 



and 



Ft 



fw.r 



£ 



fw 



( 3 - 

V 2 



dr 



sin 



fw . 2 , 
Sill I 



Uf 



(17) 



(18) 



The term 7d in eq. (0} is defined as the damping rate with 
respect to normalized amplitude, Sw = Sv/vf. Substitutions 
of eqs.l^) and ( 1181 into eq. ( 1161 gives an expression for the 
variation Sw in the stellar winds. 



dr 



{Sw) 



-Sw 



Id 



1 



1 dp 1 dvf 
v r + Vf 2p dr Vi dr 

' (3 - sin 2 6)v r + 2u f 



2{v r + Vf) 



d_ 

dr 



(- + sin 2 6)v r + Vf 



(19) 



The first term on the right hand side denotes the collisionless 
damping. The second term indicates the amplification that 
occurs in a stratified atmosphere in which the density p is de- 
creasing. The third term appears because Vf is used for the 
normalization of Sw. The fourth term is due to geometrical 
expansion of flow tubes and the variation of Vi 

For a given wave frequency, ui, the damping rate 
7d (= 7d(/9) is determined by the plasma j3 (from eq. 

|2ji and the angle 6 between the propagation direction and the 
field line, for the assumed radial propagation, that is defined 
by 

Ba 



= tan 



B , 



(20) 



Thus we see that the damping rate, 7<j, is fully determined by 
the background wind properties. 

The initial conditions of the outgoing fast mode waves, at 
the surface of the star, provide the amplitude, Swq = Svq /v f j o, 
and the period, r(= 2n/w). In the following calculations, we 
consider the results for the specific initial amplitude Swq — 
0.1. In a region close to the surface, the plasma is generally 
collisional, i.e. A(= VfT) < l m f p . In this region, eq. dl9> , for 
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collisionless damping, is not applicable, so we assume 7d = 
when solving eq. (1191 . However, even without the damping, 
Sw decreases because of the rapid increase of Vf , (see f]4). As 
a result, all other damping processes should not be effective 
either, so our treatment has acceptably isolated the collision- 
less damping case. Once the collisionless condition, eq. (|5}, 
is satisfied at r = r c \ and beyond, we can start the outward 
integration of eq. dl9l to determine Sw. Using this, we can 
estimate heating, Q w , and find the acceleration, —^§jr, that is 
caused by the dissipation of the fast MHD waves. 

4. RESULTS 

For our model of collisionless dissipation of fast MHD 
waves in a magneto-rotational wind, there are two impor- 
tant requirements. First, a sufficiently large rotation rate is 
required to create a sufficiently large 9 for the wave dissi- 
pation to occur. Second, stellar winds with moderately low 
plasma (3 values, (i.e. (3 ~ 0.01 — 0.1) are favored, because 
the fast waves, which are essentially magnetic waves, can be 
dominant only for low (3 plasmas, although the dissipation is 
suppressed if the (3 value is too small, (eq. 0}. 

4.1. Solar-type Stars 

Moderately low (3 winds are found in the coronal winds of 
low to intermediate mass stars. Although the Sun is certainly 
one of the candidates in the sense that it is a moderate to late 
type star, it has a rotation rate that is too slow to produce the 
large 9 assumed in our discussion. Observations of F-K (mid- 
to late-type) stars in open clusters (lBarnesll2003l) show that 
stars with age < 100 Myr have as large spin 5 as So — 0.5. 
Thus these are good classes of stars for us to be considering 
here. 

Figure 0] presents result of moderate magneto-rotational 
wind for So = 0.15 in a solar-mass star (M = 1M Q and R ~ 
li?o). We adopt B rfi = 8G and T = 2 x 10 6 K as typical val- 
ues. A larger mass loss rate, Me = 2 X lO _13 M0yr , is em- 
ployed than the present solar value, M ~ 2 X 10 _14 M Q yr _1 , 
because young solar-mass stars ar e found to have larger mass 
loss rate (IWood et all2002[ 2005) possibly as a result of their 
enhanced coronal activities or stellar rotation. 

On the left of fig|4] we show the wind structures as a func- 
tion of r. The top panel shows the velocity distributions, v r 
and Vff, of our chosen empirical model, and compares them 
with both vj and the rigid rotation speed rfi. The Michel ve- 
locity is vm = 1018 km/s, and is selected as following eq. 
J10I . This value is larger than the escape speed, u OS c,o = 617 
km/s. For the radial velocity distribution, we adopt the veloc- 
ity law index r\ = 3 (in eq[7J, although we find that r\ only 
weakly affects the results. It is the difference between v$ and 
the solid body rotation, rfi, leads to the spiral pattern. The 
middle panel shows that the angle 9 increases monotonically 
up to > 70°. The bottom panel presents the plasma (3 value 
(eq. 10), and it shows that moderately low-/3 circumstances 
are in fact realized for these calculations. 

Once we have chosen the empirical wind structure for a star, 
we can solve for the propagation of the fast waves by using 
eq. \\9\ . On the right side of fig. 0] results for wave peri- 
ods t = 1, 10, 100s. The top panel exhibits variation of Sw 
(solid) in comparison with the no-dissipation case (dashed). 
Although Sw decreases initially owing to the rapid increase 
of Vf to satisfy the wave energy equation (I16> . farther out it is 

5 S = 0.005 for the sun. 
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FIG. 4. — Wave dissipation in a solar-type star with So = 0.15. All the 
physical quantities are as a function of stellar radius. Upper left : velocity 
structure. Middle left : 8. Lower left : plasma f) value. Upper right : Nor- 
malized amplitude of the fast waves. Solid lines are results with dissipation 
for three cases of wave period, r. Dashed line is a result without dissipation. 
Middle right : Comparison of heati ng, Q w (erg g _1 s~ 1 ), for the three r's 
(solid) with adiabatic cooling in eq. 1211 (dashed). Lower right : Comparison 
of r. Solid lines are wave acceleration, T w .dashed line is magneto-rotational 
acceleration, T^, and dot-dashed line is acceleration by gas pressure, F g . 



amplified because of the outward propagation through a de- 
creasing density region. Waves with shorter periods, r, suffer 
greater damping per unit length, so that they dissipate nearer 
the star. 

The middle-right panel of fig|3] shows a role of the wave 
dissipation on the wind energetics. An equation describing 
variation of internal energy per unit mass, e, is written as 



de 
dr 



dr 



(r 2 v r ) - Q w = 0, 



(21) 



where the second term denotes the adiabatic loss rate and the 
third term is the heating caused by the wave dissipation. In 
our current study we evaluate the second and third terms using 
the chosen empirical wind model. We can then determine the 
magnitude of the heating term for different r (solid) and com- 
pare these with the adiabatic cooling rate (dashed). The panel 
illustrates that the region heated by the waves depends on the 
value of t. Waves with the smaller values of r heat the inner 
regions. Note that the heating exceeds the adiabatic cooling 
in the respective regions, and this indicates that the assumed 
temperature (= 2 x 10 6 K) could indeed be maintained by the 
wave dissipation. In the case of actual stars, the heating could 
be distributed more uniformly than shown in the figure, if we 
were to take into account a spectrum of waves an d inc lude 
thermal conduction as well, as will be discussed in ( i|5.1> . 
The lower-left panel of fig|4]shows the wind dynamics. The 
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equation of momentum density is 



S(deg) 
90 



dr [ 2 



1 dp g 1 dp fw GM+ 
p dr p dr r 2 



dr 4npv r 

(22) 

( Lamers & CassinellilH999i) . and the second and third terms 
are accelerations by the gas and the wave pressure, respec- 
tively. A fifth term accounts for magneto-rotational accelera- 
tion and it consists of centrifugal force and acceleration owing 
to the slinging effect of the rotating magnetic field. We com- 
pare these three terms in ratios relative to the gravitational 
acceleration, i.e. 



1 

p dr 



GM* 



1 dpfw GM* 
p dr r 2 



dr AnpVr r 2 



The magnitude of these terms has been evaluated using the 
empirical wind model and the results are shown in the panel. 
Note that the gas pressure and the magneto-rotational force 
are almost comparable in the outer region, r > lOi?*. This 
shows that it is reasonable to call the wind "a moderate 
magneto-rotator". The figure also illustrates that the wave 
pressure becomes as effective an acceleration term as these 
other two components for a magneto-rotational wind. This 
means that the fast wave damping is equally important in the 
wind dynamics; thus the wave pressure could further acceler- 
ate the wind. 

Next, we study effects of rotation rate. In fig. [5] we 
compare cases of slow rotator (So = 0.01) and fast rotator 
(So = 0.5). The angle 9 is greatly different in both cases as 
illustrated in a left panel. Accordingly, the waves dissipate 
more close to the star for the fast rotator (a middle panel), 
and the heating (and acceleration not shown) for the fast ro- 
tator occurs in the inner region (a right panel). As a result, 
the wave energy effectively dissipates in the case of the fast 
rotator to heat and accelerate the plasma, while for the slow 
rotator, the contribution is not as large. This indicates that 
regions heated by the wave dissipation depend on the stellar 
rotation even when the same waves amplitudes are incident at 
the base. 

As a star evolves, the spin rate So generally decreases. Us- 
ing this tendency, we can speculate on how stellar evolution 
affects the wave dissipation. Observationally, it is known that 
young solar-types stars have both a more rapid rotation and 
show a stronger X-ray luminosity than is the case for the sun 
iPizzolato et al.120031) . This indicates that some processes are 
operating to heat the coronae in the younger, more rapidly ro- 
tating stars, in addition to the "basal" heating processes occur- 
ring in the present solar corona. We suggest that the process of 
collisionless damping of fast MHD waves can play an impor- 
tant role in the enhanced heating of the coronae and winds and 
the X-ray source regions of young solar-type stars because it 
is much more effective in the rapidly rotating young stars. 

4.2. Massive Stars 

We have shown that our process possibly works in low 
mass main sequence stars with moderately fast rotation. In 
this subsection, we would like to consider cases of massive 
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FIG. 5. — Comparison of the propagation angle 6, amplitude 8w, and heat- 
ing rate Q w (erg g —1 s" 1 ) for a fast and slow rotating solar-mass star. The 
solid lines show results for the fast rotator, So = 0.5, and the dashed lines 
refer to the slow rotator, So = 0.01. Note that the angle and heating rates for 
the fast rotator are larger, but the wave amplitude (middle panel) for the fast 
rotator is smaller. This is because it is the radial decrease in the amplitude 
that leads to the larger heating shown in the right panel. 



stars. Traditionally, massive main sequence stars have not 
been thought to have activities concerning magnetic fields, 
i.e. coronae and/or flares, because their interiors have an 
outer radiative zone which impedes the rise of magnetic fields 
and there was no known dynamo mechanism that would 
lead to surface fields. Nevertheless, as we have stated in 
^ there has been increasing observational evidence for sur- 
face fields on hot stars. Coronal activities have been con- 
sidered to explain the X-r ay properties of mid to late B-type 
near-m ain-sequence st ars JCo hen. Cassinelli. & MacFarlane I 
119971) . More recently. iMacGregor & Cassinelh N2003I) have 
shown that the massive stars can have a surface field as a result 
of a buoyant rise through the radiative zone of fields generated 
at core-envelope interface dynamo. Although the coronal like 
geometry might be rather different from the less massive stars, 
we consider it interesting and instructive to apply the same ba- 
sic model to massive stars. 

It is known from observations of spectral lines that 
early type stars are commonly fast rotators. This is be- 
cause such a star tends to be rotating quickly after its ini- 
tial collapse and formation, the lifetimes are rather short, 
and thus the massive stars are able to persist in hav- 
ing rapid rotation. The observed mass-loss rate, M bs, 
of these types of stars is known for given M+ and 
( de la ger. Nieuwenhuiizen. & van der Huchtl[T98l . Gener- 
ally, Me ^ -Mobs in the magneto-rotational winds since equa- 
torial mass flux is larger than polar mass flux. Magnetic activ- 
ities and/or stellar rotation might further enhance intrinsic M 
as discussed in the previous subsection. To take into account 
these effects, we consider Me — (1 — 10)M o b s - Then, we 
can construct the w ind structure for a certain star after fixing 
B r and So ( §3. II . Here, we would like to introd uce a char- 
acteristic mass-loss density JCassinelh et all2002l) defined as 



M E 



(23) 



We have found that p c is a good indicator for our process; the 
damping of the fast waves work similarly in winds with the 
same p c but different M*, i?*, and Me- This is because the 
damping is controlled by density in the winds by way of the 
plasma (3 value which appears in eq.® for the damping rate. 

We now examine the roles of the fast waves in the wind en- 
ergetics and dynamics for three stellar parameters, p c , B r ,o 
and So- Since our process does not work dominantly in 
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FIG . 6 . — Critical magnetic field for collisionless damping of the fast waves 
as a function of characteristic density, p c for spin parameter So = 0.15 
(thick) & 0.5 (thin). We here consider the waves with r = l(s) and 5wq = 
0.1. Solid lines indicate criteria based on the energetics and dashed lines are 
those of the dynamics. Upper characters on the top denote corresponding 
main sequence stars (see text). In upper left region on each line the wave 
dissipation gives larger contributions to the heating or acceleration than the 
other processes. In upper right region which is indicated as 'Collisional', the 
plasma is still collisional even at r = 2R* due to high density so that our 
process does not work effectively. 

the slow rotators because the value of 9 is too small in the 
winds, we study cases with So = 0.15 (moderate rotator) 
and 0.5 (fast rotator). We investigate the wave propagation 
in main sequence stars with M* = 1 — 8M©. R± for a given 
M± is determined by a relatio n for the main sequence stars 
llKippenhahn & Weigert lU990l) . Temperature is assumed to 
be T = 2 x 10 6 K, independent of M*. Several cases of 
Me for a fixed (M*, i?*) are calculated within the range of 
Me = (1 — 10)A/ O bs- We consider waves with r = Is and 
Swq = 0.1 and investigate how the wave dissipation affects 
the energetics and dynamics for different B r fl. 

As for the energetics, we examine whether the heating in- 
tegrated with respect to a certain region exceeds the adiabatic 
cooling integrated in the same region, namely, we determine 
the condition of B r0 which satisfies 

rr ou t /-rout ^ 

/ drpr 2 Q w > / drpa 2 — (r 2 v r ). (24) 

Jr. Jr. dr 

We adopt r out = lOi?* to focus on the heating in the inner 
wind region. We also carry out the similar procedure for the 
dynamics. We compare the acceleration due to the wave pres- 
sure integrated in the same region with the acceleration by the 
other two components: 

— / drv r r 2 — > — / drv r r 2 -^- 
Jr, dr Jr. dr 

/■ r «* 2 d rQBrfy .... 

+ / drpv r r — (— ^) (25) 

Jr, dr 4TTpv r 

Figure [6] shows the critical values for B r o with respect 
to the energetics (solid; eq. ( 1241 ) and dynamics (dashed; 



eq. i25\ ) as a function of p c in the moderate (thick) and 
fast (thin) rotators. The collisionless damping of the fast 
modes can dominate the other processes if the field strength 
exceeds the respective lines. If M b s is used when deriving 
p c instead of Me in eq. d23i . p c has one-to-one correspon- 
dence to M+ because dependence of M n bs o n M+ is steep 
dde Jager. Nieuwenhuiizen. & van d er Hucht| 119881) while 
those of R+ and are weaker llKippenhahn & Wei gert I 
119901) . At the top of the fi gure, we give the spectral type of 
the stars used to determine p c . Please note that this is based 
on Mobs- If a star has larger Me as discussed previously, its 
p c becomes larger than the shown spectral type. 

Figure [6] shows that the required value for B r $ is < 100G 
even in mid to late B stars. This is a reasonable range for these 
types of stars as has been shown by Maheswaran & Cassinelli 
(1992). Cassinelli et al.(2002) have shown that such fields 
are sufficient to play a role in the channelling and torquing 
of the moderately weak winds of B main-sequence at about 
spectral type B3V and later. The critical B r o is a monotoni- 
cally increasing function of p c and a relation, B 2 /p c ^const, 
holds. Since we assume the same gas temperature, the rela- 
tion can be understood that the similar plasma (3 is required 
even for different stellar mass. As p c increases, the collision- 
less condition breaks down in the inner regions. On the dot- 
ted lines the plasma is still collisional at r = 2i?* for the 
waves with r = Is. From the magnetic fields, of order 100 
Gauss or so, that are required, we can conclude that our pro- 
cess can operate effectively in stars with spectral type of about 
B3V and later with hot outer atmospheres/winds. Interest- 
in gly, this spectral type is consistent with th e result, obtained 
bv lCohen. Cassinelli. & MacFarlane I ( 119971) . that the X-rays 
from mid- to late- B type stars seem to be arising increasingly 
from the presence of extensive hot winds and confined coronal 
zones. 

In both cases of So = 0.15 and 0.5, the condition that the 
damping dominate the dynamics demands a larger B r ,o than 
is the case for the energetics. This is because the dynamics 
condition compares the wave acceleration with both thermal 
and Poynting processes while the energetics condition com- 
pares the damping with the thermal processes (i.e. adiabatic 
cooling) only. In a region between the two conditions, the 
wind could be heated by the fast waves, but still accelerated 
dominantly by the magneto-rotational force rather than by the 
waves. For the moderate rotators, the wave dissipation is 
slower and a certain fraction of the wave energy remains at 
the radius r out . Therefore, the magnetic field B r .Q that is re- 
quired for the energetics condition is larger than that for the 
fast rotators, and this is to increase the initial wave energy. 
On the other hand, the dynamics condition shows that the fast 
rotators demand the larger B r ,o, and this is to let the wave ac- 
celeration exceed the magneto-rotational acceleration, which 
is also larger for the fast rotators. 

5. DISCUSSION 

So far we have studied outward propagation of the fast 
waves based on several assumptions : (i) wave amplitude 
at the surface (6wo = 0.1) (ii) monochromatic wave spec- 
trum (iii) no refraction (iv) a fixed, pre-specified magneto- 
rotational wind. First, we discuss these assumptions. Then, 
we compare the collisionless damping with the steepening 
damping and with other types of waves as well. Finally we 
mention another application of our process, to pulsar wind 
nebulae. 
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5.1. Limitations of Model 

(i) Wave Generation : In fj^we have assumed that Swq = 0.1 
at the inner boundary without specifying processes of the 
wave generation. The inner boundary coincides with not the 
photosphere but the base of the corona, since we assume 
T = 2 x 10 6 throughout the atmosphere. The density gen- 
erally decreases by several orders of magnitude from the pho- 
tosphere to the corona, and the region between them is called 
chromosphere and trangition region. 

We can consider two types of the wave generators. One type 
of wave generator is associated with the granulations originat- 
ing from the surface convection zones which exist in the low- 
to intermediate stars. The other type of generator corresponds 
to oscillations of closed magnetic loops emerging from the 
surface. The up-going waves generated by the former pro- 
cess suffer damping in the chromosphere before reaching the 
corona (see Stein Schwartz 1972 for sound waves). The 
closed magnetic loop process directly excites the waves at 
higher altitude in the corona (ISturrocklll999i) . and thus is 
more suitable for the present paper. Stars with faster rota- 
tion are expecte d to have stronger activities involving these 
magnetic loops dJardine et al.l2002l) . 

Accordingly, Swo would be larger in these stars than in stars 
with slower rotation; probably Swq = 0.1 is realizable in 
the fast and moderate rotators by the direct excitation in the 
corona, while it is less likely in the slow rotators due to in- 
sufficient activity by the loops. In addition to the difference 
of the dissipation conditions (fj4}, the difference of the wave 
generation processes would be important when discussing the 
stronger X-ray activities in faster rotating stars. 

(ii) Wave Spectrum : Although we have only considered the 
monochromatic waves in the model calculations, they are 
likely to be injected from actual stars with some spectrum 
of frequencies, and the wave frequency controls the wave 
damping as we have shown in the previous section. Namely 
waves with higher frequencies dissipate in the inner region. 
If some spectrum (e.g. a power-law in lo) were given for the 
input waves, the stellar atmosphere would be heated more uni- 
formly; with the inner regions is heated by waves with higher 
frequencies, and outer regions is heated by those with lower 
frequencies. 

(iii) Wave Refraction : We have assumed that the fast mode 
propagates radially without being refracted. In fact, the re- 
fraction in not negligible unless the plasma (3 is sufficiently 
small. This is because the anisotropy of Vf appears in cases 
in which the sound speed is not sufficiently small relative to 
the Alfven speed (eq. (I14» . Generally, a wave is refracted 
toward a region with a smaller phase speed. For a fast mode 
wave, it is refracted in a direction away from the field line (to- 
ward larger 9). This effect causes the fast wave to dissipate 
more quickly since the damping increases with the angle 9. 
As a result, we should expect the heated region to be shifted 
slightly inward as compared with the heated region from the 
radial propagation that we assumed. 

(iv) Fixed Background : In this paper, we treat the wave propa- 
gation within a given empirical wind structure. We think that 
this approach is acceptable as a first attempt to assess the im- 
portance of heating and acceleration produced by wave dis- 
sipation. Needless to say, a self-consistent treatment of the 
wind structure would be required in a more complete mod- 
eling of a star. This is because wave dissipation can directly 
affect both v r (via acceleration) and T (via heating). These 
two variables affect other quantities associated with the out- 



flow as well. However, we are of the opinion that empirical v r 
that we have assumed gives reasonable wind structure, simply 
because the velocity distribution is determined mainly by the 
basic stellar parameters M*, i?*, Sq. Figure[6]also illustrated 
that a stronger field was needed to modify the dynamics than 
is the case for the energetics. 

Given the velocity structure, it is the energetics and the re- 
sultant temperature that must be considered more carefully. 
Unlike the wind velocity, the temperature structure is not de- 
termined by the basic stellar parameters but rather by the 
boundary conditions for the wave generation and the subse- 
quent wave heating and radiative cooling and thermal conduc- 
tion cooling . An important point for our collisionless process 
is that the damping rate is controlled by temperature through 
(3 value (eq. (0}). Once the plasma is heated, i.e. (3 increased, 
by the waves, the damping is enhanced. It leads to further 
heating until most of the wave energy dissipates. This seem- 
ingly catastrophic behavior is partially inhibited in reality by 
the cooling effects of thermal conduction which become more 
effective if the temperature increases. In more intensive stud- 
ies, it would be important to examine the more complicated 
energetics with the other effects are accounted for. 

5.2. Collisionless Damping vs. Steepening 

We would like to compare the two types of the damping dis- 
cussed earlier, in light of the magneto-rotational wind results. 
Figure |3] shows that the collisionless process and steepening 
work in different regimes, namely the former is important for 
larger 9 except 9 = 90° and smaller Sw < 0.1 while the latter 
dominantly works for smaller 9 and larger Sw > 0.1. 

In a stellar atmosphere the amplitude, Sv, of the upgoing 
wave is inevitably amplified because the density decreasing. 
Even the normalized amplitude, Sw, increases after the slight 
decrease in the inner region as shown in the top right panel 
of Figure 0] In the cases of moderate and fast rotators, 9 is 
already large when Sw becomes large enough for the dissi- 
pation. Therefore, the collisionless damping sets in first. Sw 
is kept small by the collisionless damping in the entire re- 
gion, which is not suitable for the steepening. In the case 
of slow rotators, the collisionless damping is not so effective 
even though Sw increases to the same level, because 9 is still 
small. As a result, the dissipation of waves through shock 
steepening is expected to be more important in slow rotating 
stars. 

5.3. Other Modes 

In this paper we have focused on the fast mode waves. 
In reality, the three MHD modes couple and have a more- 
or-less turbulent nature. Since the waves and turbulence 
are mainly excited from the stellar surface, one has to con- 
sider the imbalanced cascade of tu rbulence in such a situation 
JCho. Lazarian. & Vishniadl2002l) . Also, because of the den- 
sity stratification, the amplitude of the turbulence would not 
decrease with radius even in the outer region, far from the 
surface. 

The slow wave has not been considered as a dominant 
heating source in the solar corona, because the slow waves 
that are generated from t he photosphere dissipate too rapidly 
dStein & Schwartzlll972l) . However, it has been shown that 
strong driving of the slow mode could occur at higher alti- 
tude as a result of the dynamical motion of the closed loops 
JSturrock|[T999i) in young sol ar type-stars that exhibit strong 
stellar activities jAvres 120001) . I n such a situation, supersonic 
slow modes can be important jBeresnvak. Lazarian. & Choi 
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120051) . Slow mode waves suffer TTD, similar to that affecting 
fast mode waves, and thus contribute to the heating. This will 
be studied in a fo rthcoming pap er. The steepening would be 
important as well dSuzukil2002l) . 

The Alfven wave is also a candidate for the heating in 
the low-/3 solar and stellar plasmas. The Alfven waves 
dissipate in densi ty-stratified atmosphere b y processes such 
as phase mi xing (iHey vaerts & Priest 119831). nonlinear mode 
conversion dCho & Lazariad 120021 iCho & Lazarian I 12001 
ISuzuki & Imrtsuka 2005). turbulent cascade triggered by re- 
flection dOughton et al.l200ll). High frequency Alfven waves 
suffer ion-cyclotron damping ( Cranmer. Field. & Kohl 1 19991) 
as well. The dissipative nature of the Alfven waves, which 
propagate along the field lines in the magneto-rotational 
winds is less modified for cases in which the wave length is 
shorter than the curvature scale of the winding field lines. So, 
Alfven waves can operate in both slow and fast rotators and 
also contribute to the basal heating in stellar coronae. The pro- 
cesses directly linking to the stellar rotation as shown in this 
paper would give additional contributions to lead to stronger 
coronal activities in rapidly rotating stars. 

5.4. Application to Pulsar Winds 

So far we have focused on the collisionless damping of the 
fast MHD waves in normal stars. This mechanism can work in 
any situation in which the fast waves travel in a low (3 plasma 
with curved magnetic fields. A pulsar wind is one of the can- 
didates. A pulsar magneto-sphere is a low (3 and collisionless 
plasma, and the MHD approximation is acceptable because 
the particles can interact with the waves through the magnetic 
fi elds. 

iLvubarskvl J2003I) considered dissipation of the fast waves 
propagating radially in pulsar winds by multiple shocks as a 
result of the steepening in order to investigate dissipation of 
the electromagnetic energy. Because the toroidal magnetic 
field component dominates owing to the rotation of the pul- 
sars, he studied the fast waves traveling perpendicular to the 
field lines, i.e. 6 = 90°. For this specific angle, the colli- 
sionless damping is not efficient, as shown in SJ2] so that the 
steepening would in fact be the more effective process in the 
wave dissipation. However, if 9 is not strictly = 90°, say be- 
cause of perturbed magnetic fields or other effects, collision- 
less damping would again become important. Although be- 
yond the scope of the present paper, relativistic effects should 
be taken into account in a proper study of collisionless damp- 



ing in pulsar winds. 

6. SUMMARY 

For the case of magneto-rotationally driven winds, we have 
considered the heating and acceleration of a stellar plasma 
caused by the collisionless damping of fast MHD waves. Fast 
mode waves propagate almost radially even in spiral magnetic 
fields thanks to the isotropic character of the propagation, 
provided that the (3 is low. The plasmas are typically col- 
lisionless, so the waves undergo collisionless damping. We 
find that collisionless damping of fast MHD waves is effec- 
tive when (i) the angle between the wave propagation and the 
field line is large and (ii) the plasma (3 is around 0.01 — 0.1. In 
the magneto-rotational winds, the angle increases as the wave 
travels away from the surface. This angle change eventually 
leads to a dissipation of the waves, and causes heating and 
acceleration of the surrounding medium. 

We have found that our collisionless damping is impor- 
tant in moderate and fast magneto-rotators of solar-mass stars, 
while it is not so effective in slow rotators such as the (present) 
sun. It cooperates with the "basal" processes which are re- 
sponsible for the heating of the corona in the slow rotators. 
We propose that the mechanism of collisionless damping of 
fast waves can explain the observed strong coronal activities 
in young (age < lOOMyr) solar-type stars because of their 
rapid rotation. The damping also effectively applies in the 
energetics, and perhaps also to the dynamics, of the atmo- 
spheres of moderately massive main sequence stars of which 
are of spectral types about B3V or later and which have ob- 
servationally plausible magnetic fields (< 100G). 

When compared with the shock dissipation by the steepen- 
ing, the collisionless damping is expected to be more impor- 
tant in moderate and fast rotators with an input wave ampli- 
tude, Swo < 0.1. In addition, the process is also argued to 
be important in general astrophysical low (3 plasma which has 
curved magnetic fields, such as pulsar magneto-spheres. 
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APPENDIX 

DAMPING BY STEEPENING 

In this appendix, we summarize the damping of fast shock trains which form by the steepening of fast MHD waves. First, 
we derive entropy generation at the fast shock whose geometry is shown in fig. Physical quantities in the upstream region 
are indicated by subscript "1" and those in the downstream region by "2". Components parallel with shock normal are denoted 
by subscript "||" and those perpendicular to it by "_L". Because of the conservation of magnetic flux, V • B = 0, the parallel 
component of magnetic field, Bn, is constant. The perpendicular field strength is increased from B± to B± + 5B± at the fast 
shock. The equations for the conservation of mass flux, momentum parallel and perpendicular to shock normal, electric field, and 
energy are expressed in the frame of an observer co-moving with the shock as (e.g. Priest & Forbes 2000) 

PlV\\,l = P2V\\,2, (Al) 

^2 ^2 , 2B ± 6B ± +6Bj 

I! II ^ 



B\\SB ± 

47T 

v\\ il B ± -v± tl B ll =VH 2 (B ± +8B±)-v±flBJ=0) (A4) 



Pl%l^_L,i = P2Vn 2 VA_.2 -. , (A3) 

" " 4-7T 
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and 



B l B \\ B ± 1 / 2 2 ^ 

PlV\\,i + Pieit'||,i + ^T u ll,i i^~ Vj -' 1 + 2 />1 ^H' 1 + ^-L.i^H - 1 



(B^+SB^) 2 B {l (B ± +5B ± ) l 

= 2 + P2e 2 W||.2 H "j W||,2 : V± 2 + — P2C^|| 2 + W _L 2) U ||,2- (A5) 

47T 47T 2 

Here we define all the quantities in a frame co-moving with the shoc k fron t. Then, both left- and right-hand sides of ea. JA4> 
become 0. By eliminating uim , vii 2 , and j<2 from eqs. JAl> - (IA5> . we have a relation between SB± and compression 

ratio, 

cj= P —, (A6) 
Pi 

as 

_ _ ^ + + 1 _ _ SB^ 1)B 2 - _ 

where £? 2 = By + The pressure jump is also derived from eq. (IA2> : 

Ap = P2 - Pl = ^ — ^ y x (A8) 

Entropy generation is written as iMihalas & Mihalas 119841) 

£ = ^H P ~4) - "^{Ml + ^) -7lnWh (A9) 
i? 7-I pip 2 ' 7-1 p 

whe re is a gas c onstant, and we assume that the shock is weak. Then, As is derived from a given amplitude, SB±, by using 
eas.llA"7l and dA8l . 

The steepening of the wave fronts forms the fast shock train. The wave dissipates by these multiple shocks. The heating rate, 
Q s (crg cm'V 1 ), of the fast shock train with angular frequency, ui, is given by 

„ As „„ uj As lu 

Q » = ll^ RT ^ = ll^ (A10) 

SB 2 

The damping rate is derived from Q s and energy density of the wave, -s^, as 



a 

6B 2 ± /8ir 

Note that if 5B± / B is fixed for the weak shocks, 7 S is independent of plasma (3 because 



7s = T7^Z (AH) 



As 1 Ap 1 2 1 

-R Pl 5Bj^^ Pl 5Bj°< B 5Bi 



7s « -frPlTT^T - —Pi T^T « B TnT ( A12 ) 



from eqs.( IA8> — JA1 It 
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